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Abstract 

In warped supersymmetric grand unification, XY gauge particles appear near the TeV 
scale along with Kaluza-Klein towers of the standard model gauge fields. In spite of this 
exotic low-energy physics, MSSM gauge coupling unification is preserved and proton decay 
is naturally suppressed. In this paper we study in detail the low-lying mass spectrum 
of superparticles and GUT particles in this theory, taking supersymmetry breaking to be 
localized to the TeV brane. The masses of the MSSM particles, Kaluza-Klein modes, 
and XY states are all determined by two parameters, one which fixes the strength of 
the supersymmetry breaking and the other which sets the scale of the infrared brane. A 
particularly interesting result is that for relatively strong supersymmetry breaking, the 
XY gauginos and the lowest Kaluza-Klein excitations of the MSSM gauginos may both 
lie within reach of the LHC, providing the possibility that the underlying unified gauge 
symmetry and the enhanced N = 2 supersymmetry of the theory will both be revealed. 



1 Introduction 



One of the most striking results of any extension of the standard model is the unification of gauge 
couplings in theories with a supersymmetric desert above the TeV scale. Certain mysteries of the 
standard model, such as the stability of the Higgs potential and hypercharge quantization, can 
be elegantly addressed by a combination of low-energy supersymmetry and grand unification at 
high energies [1], making the supersymmetric desert seem even more compelling. An immediate 
consequence of this framework is the presence of superparticles at the electroweak scale. In the 
standard paradigm these superparticles are assumed to be in the smallest possible supersym- 
metric representations: the standard model particles and their superpartners form multiplets of 
N — 1 supersymmetry, and there are essentially no other fields at the TeV scale charged under 
the standard model gauge group. This is the basis of the minimal supersymmetric standard 
model (MSSM), which has been the main focus of phenomenological studies in physics beyond 
the standard model. 

In a previous paper [2] we studied an alternative to this framework with Goldberger. We 
constructed a realistic theory of grand unification in warped space, in which the unified gauge 
symmetry is broken by boundary conditions and the electroweak scale is generated by the warp 
factor. The theory predicts a rich spectrum of new particles at the TeV scale; in addition to the 
usual superpartners, there are Kaluza-Klein (KK) towers for the standard model fields as well as 
their supersymmetric and SU (5) partners. The appearance of these particles allows the theory 
to be "higher dimensional" at the TeV scale. In particular, radiative corrections to the Higgs 
potential are extremely soft, namely exponentially shut off above the TeV scale, and we can have 
a complete understanding of the MSSM Higgs sector through the U(1)r symmetric structure 
of the theory. We showed that the theory preserves the successful MSSM prediction for gauge 
coupling unification, despite the drastic departure of the matter content from that of the MSSM 
at the TeV scale. The theory also preserves a number of the successes of conventional high- 
scale unification: proton decay is sufficiently suppressed and small neutrino masses are naturally 
obtained through the seesaw mechanism. This theory thus naturally synthesizes two dominant 
approaches to physics beyond the standard model: Planck-cutoff and TeV-cutoff paradigms. 
Some of these features were suggested earlier in [3], and an understanding of logarithmic gauge 
coupling evolution in warped space was developed in [3-9]. 

In this paper we study the phenomenology of the warped supersymmetric grand unified 
theory (GUT) described above. Due to the rich spectrum near the TeV scale, the experimental 
implications of the theory can be quite different from the conventional MSSM. In the dual 4D 
picture of our theory supersymmetry is broken at the TeV scale by strong dynamics. In the 5D 
picture, the effects of this supersymmetry breaking are parameterized by operators localized on 
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the TeV brane. The phenomenology of the theory then crucially depends on the form and size of 
these operators and the location of matter and Higgs fields in the fifth dimension. In this paper 
we consider the case where matter fields are localized to the Planck brane, which is consistent 
with the requirements from proton decay suppression and gauge coupling unification. The Higgs 
fields can either be localized to the Planck brane or propagate in the bulk, though we focus on 
the localized case in the latter part of our analysis. The supersymmetry breaking operator on 
the TeV brane is taken to be a linear term for a singlet superfield in the superpotential, which 
was introduced in [10] and considered in the unified theory context in [2, 11]. This gives gaugino 
masses at tree level through an operator localized on the TeV brane. At one-loop level, squarks 
and sleptons obtain masses that are insensitive to physics above the TeV scale. These masses 
are flavor universal, so that the supersymmetric flavor problem is naturally solved in this setup. 

An important feature of the present framework is that the masses for the electroweak-scale 
particles are determined in terms of only a few free parameters. In [2] we have shown that 
the simplest theory of warped GUTs is obtained with the TeV brane respecting the full SU(5) 
symmetry and the Planck brane respecting only the standard model gauge symmetry. This 
implies that the operators on the TeV brane, including the one that generates the gaugino 
masses, must respect SU(5). The coefficient of this gaugino mass operator determines the 
masses for all the superparticles, GUT particles and KK towers, up to the overall mass scale 
and small effects from electroweak symmetry breaking. This situation is quite different from 
that in non-unified theories [10], where we can have arbitrary values for SU(3)c, SU(2) L and 
U(l)y gaugino masses and thus have less predictive power. An especially interesting parameter 
region for our theory is where the coefficient of the gaugino mass operator becomes large. In 
this parameter region, the spectrum appears quite different from that of the MSSM: the MSSM 
gauginos become pseudo-Dirac states and one of the XY gauginos becomes quite light, even 
lighter than some of the superparticles [2]. This makes future experimental searches for these 
particles quite exciting and promising. 

The organization of the paper is as follows. In section 2 we review the theory and summarize 
the framework for our computation. In section 3 we calculate the masses of the superparticles 
and GUT particles, including one-loop radiative effects. We find that even with the present 
experimental bounds on the superparticle masses, the lightest XY gaugino is well within the 
reach of the LHC for moderately large supersymmetry breaking. In section 4 we discuss the 
physical Higgs-boson mass and the naturalness of electroweak symmetry breaking. Conclusions 
are given in section 5. 
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2 Theory and Framework 



2.1 Warped supersymmetric GUTs 

We begin by reviewing the warped supersymmetric grand unified theory of Ref . [2] . The theory 
is formulated in a warped 5D spacetime with the extra dimension compactified on an S 1 /Z 2 
orbifold: < y < nR, where y represents the coordinate of the extra dimension. The metric of 
this space is given by 

ds 2 = G M Ndx M dx N = e-^ri^dx^dx" + dy 2 . (1) 

Here, k is the AdS curvature, which is taken to be somewhat (typically a factor of a few) smaller 
than the 5D Planck scale M 5 ; the 4D Planck scale, Mpi, is given by Mp! ~ Ml/k and we take 
k ~ M 5 ~ Mpi. We choose kR ~ 10 so that the TeV scale is naturally generated by the AdS 
warp factor: T = ke~ nkR ~ TeV [12]. 

We consider a supersymmetric SU (5) gauge theory on the above gravitational background. 
The bulk SU{h) symmetry is broken by boundary conditions imposed at the boundary at y — 0. 
Specifically, the 5D gauge multiplet can be decomposed into a 4D N = 1 vector superfield 
V(A fl , A) and a 4D N = 1 chiral superfield S(a + ?A 5 , A'), where both V and £ are in the adjoint 
representation of SU(5). The boundary conditions for these fields are given by 

(s) W-y) = {-pzp-i) (e) {x ">- y,) = {-x) ix "> y>) > (2) 

where y' = y — irR, and P is a 5 x 5 matrix acting on gauge space: P = diag(+, +,+,—,—). This 
reduces the gauge symmetry at the y = brane (Planck brane) to SU(3)c x SU(2) L x U(l)y 
(321), while the 5D bulk and the y = nR brane (TeV brane) respect full SU(5). After KK 
decomposition, the above boundary conditions ensure that only the 321 components of the 4D 
vector superfield, V, have massless modes. The typical mass scale for the KK towers is T ~ TeV, 
so that the lowest KK excitations of the standard model gauge fields and the lightest XY gauge 
bosons both have masses of order TeV. In fact, the KK towers for these gauge fields turn out to 
be approximately SU (5) symmetric. 

The Higgs fields are introduced in the bulk as two hypermultiplets transforming as the 
fundamental representation of SU(5). Using the notation where a hypermultiplet is represented 
by two 4D N — 1 chiral superfields $(0, ip) and <3> C (0 C , ip c ) with the opposite gauge transformation 
properties, our two Higgs hypermultiplets can be written as {H, H c } and {H, H c }, where H and 
H c transform as 5 and H and H c transform as 5 under SU(5). The boundary conditions are 
given by 

(hc) (s*,-v) = (~ p P Hc)(^,yl (|c) = ( H H c) (^,y0, (3) 
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(p,p') 


gauge and Higgs fields 


bulk matter fields 


(+,+) 
(-,-) 
(-,+) 
(+,-) 


V321, Hd, Hd 

V £7"c fjc 

^321, tl D) ii D 
Vx, H T , H t 
^x, H T , H T 


Tu,e, Tq, F d , F' l 

rpc t ,c T? c T? ,c 

Tq, T(je, F l , F' d 

rpc rplc rpc 77>/C 

1 Qi 1 U,E-> "hi " D 



Table 1: Boundary conditions for the bulk fields under the orbifold reflections. Here, Tq UE 

(Fdl) are ^e components of T^> (F^) decomposed into irreducible representations of the stan- 
dard model gauge group. The fields written in the (p,p f ) row, ip, obey the boundary conditions 
<p(-y) =P ( f(y) and ip(-y') =p'tp(y'). 

for {H,H C }, and similarly for {H,H C }. After KK decomposition, only the two Higgs doublets 
from H and H have massless modes. All the other KK modes, including those of colored Higgs 
fields, are massive with characteristic mass scale given by T ~ TeV. The masses for these modes 
are approximately SU(5) symmetric as in the case of the gauge fields. Therefore, at this stage, 
the mass spectrum of the theory is given as follows: we have a 321 vector multiplet, V321, and 
the two Higgs doublets, Hd and Hd, at the massless level together with SU(5) symmetric (and 
N = 2 supersymmetric) KK towers for the gauge and Higgs fields with characteristic mass scale 
T ~ TeV. 

Despite the fact that XY gauge and colored Higgs fields have masses of order TeV, proton 
decay can be adequately suppressed. One way to achieve this is to impose baryon number, 
which is possible even if matter propagates in the bulk, and this approach allows matter fields 
to have wavefunctions spread over the extra dimension. Another, probably more satisfactory 
way is to localize matter toward the Planck brane - either strictly localized as brane fields, or 
approximately localized using bulk mass parameters. In particular, the latter possibility appro- 
priately quantizes matter hypercharges while avoiding rapid proton decay caused by exchanges 
of TeV-scale GUT particles. 

To be more explicit, the bulk matter model has the following four hypermultiplets for 
each generation: {T, T c }(10), {T f , T /c }(10), {F, F c }(5*) and {F', F /c }(5*), where the numbers 
in parentheses represent the transformation properties of the non-conjugated fields under 577(5). 
The boundary conditions for the matter fields are given similarly to the Higgs fields, Eq. (3), 
but for {T, T c } and {T', T /c } the matrix P acts on both 577(5) fundamental indices and the 
overall parities under y — > — y are taken to be opposite between T and T' multiplets and be- 
tween F and F' multiplets (these boundary conditions are summarized in Table 1). With these 
boundary conditions, a complete generation, Q, U, D, L and E, arises at the massless level as 
T(U, E),T'(Q), F(D), F'(L). The wavefunction profiles for these modes depend on the bulk 
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hypermultiplet masses 

cttR 



s = 



J d 4 x £ dy V^G J d 2 6c^k^ c + h.c. 



(4) 



parameterized by dimensionless quantities c$, where $ runs for T,T',F and F' . For c$ > 1/2, 
we find that the wavefunctions for the zero modes are strongly localized to the Planck brane as 
e -(c*-i/2)fc|j/| anc j fasti proton decay rates are sufficiently suppressed for c$ <: 1. 

From the 5D viewpoint, there are three local operators that can contribute to the low-energy 
4D gauge couplings: 

i-kR r 1 1 11 



S =- [d 4 x r dyV^G\—F MN F MN + 25(y)—F% u F a ^ + 25(y - nR)-F, u F^ 



where the index a runs over SU (3)c, S'C/ (2) L and f/(l)y (a = 3, 2, 1, respectively). The structure 
of these terms are determined by the restricted 5D gauge symmetry, which reduces to 321 on 
the y = brane but is 577(5) at all the other points in the extra dimension. At the fundamental 
scale M* ~ M 5 , the coefficients of these operators are incalculable parameters of the effective 
field theory. Therefore, one might worry that one cannot obtain any prediction for the low- 
energy gauge couplings, which in general depend on these unknown parameters. This difficulty, 
however, is avoided by requiring that the theory is strongly coupled at the scale M*. In this case 
the sizes of these coefficients are estimated as l/g% — M*/167r 3 and \ jg\ a — l/g% — 1/IQtt 2 , and 
one finds that the low-energy prediction is insensitive to the parameters go,a and g n evaluated 
at M*, which encode unknown physics above the cutoff scale of the theory. The prediction for 
the low-energy 4D gauge couplings, g a , is then written in the form 

1 ~ (£77(5) symmetric) + ^A a (T, k), (6) 



gl(T) v v ' " 87r> 

where A a (T, k) is the quantity whose non-universal part can be unambiguously computed in the 
effective theory. In our theory, this quantity is given by 

/ u\ 

(7) 

at one-loop leading-log level, regardless of the values of the bulk mass parameters as long as 
they are larger than or equal to 1/2: ch,Cj},ct,ct>,cf,cf' > 1/2 (we have absorbed a possible 
SU(5) symmetric piece into the first term of Eq. (6)). This is exactly the relation obtained in 
conventional 4D supersymmetric unification with the parameter k identified with the unification 
scale. Therefore, we find that our theory preserves the successful 4D MSSM prediction for gauge 
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A 2 (T, k) ~ 
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coupling unification, despite the drastic departure of the matter content from the MSSM at the 
TeV scale [2]. 

Given that MSSM gauge coupling unification is naturally preserved when matter and Higgs 
fields have c > 1/2 and that proton stability is ensured when matter fields have c ^ 1, it is 
natural to focus on the case with all matter fields strongly localized to the Planck brane. The c 
parameters for the Higgs multiplets are less constrained, but we can certainly consider the case 
where the Higgs fields are also effectively localized to the Planck brane. In this case the physics 
is well approximated by simply regarding matter and Higgs as brane fields, as will be done in 
our calculations in the subsequent sections. (The analysis of section 3 assumes matter to be 
localized to the Planck brane but is independent of the Higgs profiles; parts of section 4 assume 
that the Higgs fields are also localized to the Planck brane.) However, one should keep in mind 
that our analyses also apply for bulk matter provided the lowest KK modes are localized toward 
the Planck brane by bulk hypermultiplet masses. 

Finally, we discuss the Yukawa couplings. The Yukawa couplings are written on the Planck 
brane as 

S = Jd 4 x J^dy V^G 25{y) J d 2 6 (y u QUH D + y d QDH D + y e LEH D ) + h.c. . (8) 

Since the gauge symmetry on the y = brane is only 321, we do not have unwanted 577(5) mass 
relations such as m s /md = m^/m e . The above Yukawa couplings respect a U(1)r symmetry, 
under which the 4D superfields V, E, H and H are neutral, T, T c , F, F c , T', T /c , F' and F' c have 
unit charge, and H c and H c have charge +2. This U(1)r forbids dangerous dimension four and 
five proton decay operators together with a potentially large supersymmetric mass term for the 
Higgs fields, thus providing a complete solution to the doublet-triplet splitting and proton decay 
problems (the U(1)r symmetry is broken to its Z 2 subgroup through supersymmetry breaking 
discussed in the next subsection, but without reintroducing phenomenological problems). Small 
neutrino masses can be naturally generated by introducing right-handed neutrino fields with 
the Majorana mass terms and neutrino Yukawa couplings on the Planck brane, through the 
conventional seesaw mechanism. 



2.2 Framework for the analyses 

To calculate physical quantities such as superparticle and GUT particle masses, we must specify 
how supersymmetry is broken. We also have to specify a calculational scheme for computing 
radiative effects, which are quite important for the phenomenology of the theory. 

Since any mass parameter on the y = 7iR brane of order the fundamental scale appears as a 
TeV scale parameter in the 4D picture, it is quite natural to consider supersymmetry breaking 
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on the TeV brane. Specifically, we introduce a supersymmetry breaking potential 
S = Jd A x J^dy 25{y - nR) J d 2 9d 2 BZ ] Z + (^j d 2 9A 2 z Z + h.c. 



(9) 



on the TeV brane [10]. Here, Z is a gauge singlet chiral superfield, and A z is a mass parameter 
of order the fundamental scale M* ~ M 5 . The resulting supersymmetry breaking is transmitted 
to the SU(5) sector through the following operator: 



S = [d A x [^dy ^f^G 2S(y - irR) [ d 2 9-^-ZW a W a + h.c. 
J Jo J 2M* 



(10) 



where W a represents the gauge field strength superfield for the bulk SU(5) gauge multiplet. 
(The presence of both operators in Eqs. (9, 10) breaks the U(1)r symmetry discussed in the 
previous subsection, but this breaking does not reintroduce phenomenological problems such as 
rapid proton decay [2].) When expanded into component fields, this gives the gaugino masses 
localized on the TeV brane 



S = Jd 4 x £ dyV z G25{y-nR) 



J/!x -\ a \ a + h.c. 



V ■ (ii) 

where M A = AA|/M*, and \ a is the SU(5) gaugino. After KK decomposition, this term gives 
TeV scale masses for the 321 gauginos. Since the full SU(5) symmetry is respected on the TeV 
brane, we find that all the gaugino masses are fixed by the single parameter M\ at tree level. 
The splitting among the three gaugino masses then arises through radiative effects. 

Squarks and sleptons, which are localized to the Planck brane, obtain masses at one-loop 
level. Because of the geometrical separation between supersymmetry breaking and the place 
where squarks and sleptons are located, the resulting masses are finite and calculable in the 
effective field theory. In fact, the loop integrals are cut off at the scale T and are insensitive to 
unknown UV physics. However, to find the detailed structure of the mass spectrum, for example 
that coming from the splitting of the gaugino masses, we have to include higher order effects. 

What calculational scheme should we use to compute superparticle masses including radiative 

effects? One way of computing radiative corrections in truncated AdSs is to calculate them 

directly in perturbation theory, using the KK decomposed 4D theory and retaining all KK modes 

in loops. This procedure is justified as long as the external momenta, p, are sufficiently smaller 

than the threshold for the KK towers, T [5, 8], because then the effects from unknown physics 

above the cutoff scale are suppressed by powers oip/M'^. Here is the cutoff scale on the IR 

brane, Ml = M*e _7rfcfi . However, this procedure is not quite suitable for computing radiative 

corrections to superparticle masses, as we expect to get powers of large logarithms at higher 

loop orders, (a/47r) n (ln(A;/T)) n for the n-th loop, which invalidate the perturbative expansions. 1 

1 For radiative corrections to the gauge couplings, the form of the one-loop renormalization group equations 
ensures that there are no such terms beyond one loop, if we compute corrections to 1/g 2 . 
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Although general renormalization theory relates the coefficients of these leading logarithms, in 
principle allowing their summation, this procedure requires computation of at least the lowest 
loop diagrams containing the large logarithms: one loop diagrams for the gauginos and two loops 
for the sfermions. These calculations are somewhat involved, so we do not adopt this scheme for 
computing the superparticle masses. 

As in 4D theories, leading-log effects can be taken into account using the renormalization 
group method. In truncated AdSs, the direct analog of integrating out higher momentum modes 
is to integrate out the space closest to the Planck brane. The specific procedure [13] is sum- 
marized as follows. Starting from the theory where the two branes are located at y — and 
y = ttR with the AdS curvature given by k, we can construct a theory in which the two branes 
are located at y = eR and y = 7rR by integrating out the region < y < eR. Then, rescaling all 
the mass scales of the theory as m — > m' = e~ nekR m, we obtain a theory with the AdS curvature 
given by k' (< k), which has the same IR scale, k' e -(*- € ) kR = ke~ wkR = T, and gives the same 
low-energy predictions as the original theory. By repeating this procedure, we can obtain the 
theory where k' is sufficiently close to T that there is no large logarithm in loop calculations. 
However, the theory obtained in this way contains a series of higher dimensional operators on 
the UV brane suppressed only by powers of k! <C k, whose effects on predictions are 0((T/k') n ) 
for some power n. Without knowing the coefficients of these higher dimensional operators, the 
only way of obtaining reliable predictions is to choose k' to be somewhat larger than T, but 
then we cannot really sum logarithms down to the scale T. Due to this limitation, we do not 
choose this "floating cutoff" scheme either, although one can compute superparticle masses in 
this scheme if one is satisfied with the precision in which one does not distinguish the values for 
running couplings at T and k' ~ T/e, where e sets the typical size of errors in the predictions. 

Instead of changing the cutoff, one can also sum up the large logarithms by defining the 
couplings of the theory using the sliding renormalization scale /i. Suppose we compute radiative 
corrections to the mass of a particle that is localized on the Planck brane. If we use the couplings 
appearing in the bare Lagrangian, the resulting expression contains large logarithms \n(k/T). 
However, these large logarithms can be successfully resummed if we use the couplings defined at 
the scale // ~ T, say through momentum subtraction, measured in terms of the 4D metric r]^. 
This procedure effectively corresponds to integrating out physics above the scale \x and encoding 
it into the couplings defined at //. For this procedure to work, the effective theory obviously has 
to be valid up to k, the scale below which the large logarithms are generated. This is indeed 
the case when we compute radiative corrections to Planck-brane localized quantities [5] , because 
on the Planck brane physics is essentially four dimensional up to the scale k. This implies 
that large logarithms that could invalidate perturbative expansions, i.e. ln/c/T's appearing in 
the superparticle mass calculation, are effectively resummed into the coefficients of operators 
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located on the Planck brane. Since we did not lower the cutoff scale of the theory, the mass 
scale on this brane is given by k, and the coefficients of these operators still scale by powers of 
k. Therefore, (UV insensitive) low-energy quantities, such as squark and slepton masses, can be 
reliably computed using the lowest dimension operators on the Planck brane with coefficients 
evaluated at the scale fj, ~ T. 

Calculating superparticle masses in our setup requires only the values of the brane localized 
gauge couplings, go,a and g n , at the 4D scale T. Recall that our strong coupling assumption sets 
these couplings to be ~ 4n at the scale M* measured in terms of the 5D metric (see the sentences 
above Eq. (6)). This corresponds to l/^o a(A* = — 1 /167T 2 and \jg\{p, — T) ~ 1/167T 2 in terms 
of the 4D metric, so all ln(/c/T) effects can be resummed by running down the couplings go,a 
from the scale k down to T. Since the solution to the one-loop renormalization group equation 
(RGE) takes the form 



the required couplings go, a (T) can be obtained from b a without knowing the precise values for 
the initial couplings, g 0ta (k). 

Because the non-universal part of the low-energy 4D gauge couplings (i.e. the differences of 
the three gauge couplings) comes only from the Planck-brane localized couplings g 0ta (T), we know 
that the non-universal part of b a must be the same as that of the MSSM: b a -b b = &mssm_ & mssm_ 
This allows us to write b a = b^ SSM + b, where b takes a universal value for all the 321 gauge 
groups. To determine the value of b, we focus on the U(1) Y component (a = 1). For a U(l) 
theory, the computation of Ref. [5] explicitly shows that contributions to b a are saturated by 
zero-mode fields for bulk scalars and fermions (for non-Abelian gauge fields this saturation was 
shown only for the non- universal part). We also notice that the XY gauge bosons and gauginos 
are all strongly localized to the TeV brane and thus do not contribute to b\. This is sufficient to 
conclude that 6 = 0, because the contribution to b\ comes entirely from the scalars and fermions 
in the bulk and on the Planck brane, giving bi = b^ lsSM . 

Our procedure for computing superparticle masses can now be summarized as follows. We 
first integrate out physics above the scale T (in the 4D metric) by defining the couplings at the 
sliding renormalization scale fi ~ T. Using these couplings, we compute lowest order contribu- 
tions to the superparticle masses: tree level for the gauginos and one loop for the sfermions. 
Higher loop effects are expected not to contain large logarithms because they are already in- 
cluded in the renormalized couplings, so that perturbation theory must work well. The remaining 
uncertainty arises from possible TeV-brane operators suppressed by powers of M'^ = M*e~ nkR . 
These operators, which are intrinsically incalculable in the effective theory, bring uncertainties 
of 0((T/MQ n ) in the predictions, where n is some power depending on the dimension of the 




(12) 
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operator. However, 4D Lorentz invariance implies that these corrections are at most of or- 
der (/c/M*) 2 ^ 10%, and we can trust our leading-order computations up to uncertainties of 
0(10%). 2 Keeping this remark in mind, in the next section we compute the mass spectrum of 
the theory. 

3 Spectra for Superparticles and GUT Particles 

In the supersymmetric limit, the massless sector of the theory consists of the fields of the MSSM. 
The gauge multiplets propagate in the bulk, so the MSSM gauge bosons and gauginos are 
accompanied by KK towers of massive gauge multiplets. These massive KK levels are N = 2 
supersymmetric and approximately £77(5) symmetric, with masses given at tree level by m n ~ 
(n — l/4)7r T for n — 1, 2, • • •. The lightest gaugino and gauge boson KK modes thus have masses 
m\ ~ 2.4 T in this limit. If they propagate in the bulk, Higgs and matter fields will also have 
KK excitations, but the analysis of this section will treat matter as localized to the Planck brane 
(or approximately localized with a large bulk hypermultiplet mass term). Our results will apply 
regardless of whether or not the Higgs fields propagate in the bulk, however. 

When supersymmetry is broken as described in section 2.2, the MSSM gauginos acquire 
masses at tree level and the squarks and sleptons obtain masses at one loop. Supersymmetry 
breaking also feeds into the spectrum of the KK excitations in potentially crucial fashion. In 
particular, the masses of the lightest XY gauginos and of the first KK excitations of the 321 
gauginos can be pushed well below the 2.4 T value that applies in the supersymmetric limit, 
improving the prospects for their discovery at colliders [2] . Our aim in this section is to compute 
the masses of the MSSM superparticles along with those of the lightest 321 KK states, XY gauge 
bosons, and XY gauginos, including one-loop radiative effects. 

Before getting into the details of these computations, we refer the reader to Fig. 1, where a 
schematic depiction of the effects of supersymmetry breaking on the particle spectrum is given. 
It is useful to introduce the supersymmetry-breaking parameter 

x = M\/k, (13) 

which we will take to be .01 — 10 corresponding to relatively unsuppressed supersymmetry 
breaking on the TeV brane. Once this single parameter is specified, the masses of all the particles 
in the theory are fixed up to an overall scale and small electroweak symmetry breaking effects. 
For x = (Fig. la), the massless particles are the 321 gauge bosons A 3 ? 1 , the 321 gauginos 

2 Corrections to the scalar masses induced by a supersymmetry-breaking gaugino kinetic operator on the TeV 
brane can be calculated using the gaugino propagators derived in Appendix A, by choosing appropriately. 
We find that these corrections do not significantly affect our results below. 
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(a) x = 



(b) x < 1 



m. 



3T 



IT 



m r 




Af 1 A 321 A XY 



4 21 A 321 A XY gj 



A 321 A 321 X M A XY g~J 



Figure 1: The lowest-lying masses for the 321 gauge bosons (A 321 ), 321 gauginos (A 321 ), XY gauge 
bosons (Xfj), XY gauginos (A XY ), and MSSM scalars (q, I), in the presence of no supersymmetry 
breaking (a), weak supersymmetry breaking (b), and strong supersymmetry breaking (c). Each 
bullet for A 321 and A XY represents a Majorana and Dirac degree of freedom, respectively. Dotted 
lines connect nearly degenerate SU(5) partners, and arrows indicate displacements of the mass 
levels relative to their supersymmetric positions. The states not shown here include the standard- 
model quarks and leptons, the two Higgs hypermultiplets, and the real scalar fields arising from 
the 5D gauge multiplet, which have the 321 and XY quantum numbers and are degenerate with 
the excited states of ^4 321 and X^. 
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A 321 , and the MSSM scalars q, I (and the standard-model quarks and leptons, of course). The 
first KK level, at mi ~ (3ir/A)T, includes the first A 321 KK mode and a degenerate pair of 
Majorana fermions for A 321 , which form a Dirac state. Nearly degenerate with these are their 
SU (5) partners, the lightest XY gauge bosons and a degenerate pair of Dirac XY gauginos 
A XY . When supersymmetry is broken by a small amount (a; < 1; Fig. lb), the A 321 zero modes 
pick up small (<C T) tree-level masses and the MSSM scalars pick up even smaller masses at 
loop level. Meanwhile, the degeneracy between the A 321 KK excitations is spoiled: one Majorana 
fermion's mass increases, while the other's decreases. The XY gaugino masses move in essentially 
the same way: one Dirac state becomes heavier while the other becomes lighter. Thus, as long 
as x is sufficiently small, the gaugino states are still SU (5) symmetric. In the limit of very strong 
supersymmetry breaking (x ^> 1; Fig. lc), the 321 gauginos once again form near-degenerate 
pairs, but now the formerly massless gaugino is approximately degenerate with one of the light 
KK excitations, with a mass ~ T/A. Meanwhile the other light gaugino KK excitation becomes 
nearly degenerate with one of the 321 gauginos from the second KK level, with mass ~ (5ir/A) T. 
On the other hand, the lightest XY gaugino has no nearly degenerate partner in this regime. 
Moreover, its mass approaches zero, rather than a finite limiting value, as x — > oo. Therefore, it 
is possible that this XY gaugino is quite light. 

Let us now compute the particle masses. To calculate supersymmetry-breaking effects on 
the spectrum, we adopt the framework described in the previous section. As discussed there, 
we assume that the effects of physics above the energy scale T are encoded in local operators 
that reside on the Planck brane. The operators of interest for our purposes are the Planck-brane 
localized gauge kinetic terms. We compute tree-level gaugino and gauge boson masses and one- 
loop scalar masses under the presence of these (radiatively induced) Planck-brane operators. 
This requires knowing the numerical values of the Planck-brane, TeV-brane, and bulk gauge 
couplings at the scale T. The TeV-brane coefficient 1/g 2 is assumed to be very small at that 
scale, ~ 1/16-7T 2 , and we neglect it entirely. The Planck-brane coefficients l/g 2 a are similarly 
assumed to be negligibly small at the scale k, and their values at T are obtained using the MSSM 
RGEs (see Eq. (12) and the discussion below it). Finally, the value of the bulk gauge coupling 
gB is determined from 

9 2 a (T) ~ 9 B + ~9i a (T) + ~gl{Ty { } 

where g a (T) are the 4D gauge couplings evaluated at the scale T. These are approximated by 
running the experimentally measured values at mz up to ~ 1 TeV using standard model RGEs 
and then from ~ 1 TeV to T using MSSM RGEs. Fixing in this way the numerical values of the 
various gauge couplings, we can now give results. 



12 



(i) 321 gauginos 

The equation determining the 321 gaugino masses in the presence of Planck-brane and TeV-brane 
gauge kinetic terms is presented in Appendix B. In the limit where the TeV-brane coefficient 
vanishes, the equation becomes 



J ° O) + IN Jl O) _ J o O) + 9 2 B M X J, 



(15) 



where we have suppressed the 321 index. Using the values for the gauge couplings at the scale 
T, the lowest mass solutions give the lightest gaugino masses at T, which we run down from T 
to the gaugino masses themselves using MSSM RGEs. We will be interested primarily in the 
masses of the two lightest sets of gauginos. 
(ii) MSSM scalars 

The MSSM scalars acquire masses at one-loop level due to their gauge interactions. We obtain 
these contributions to the scalar masses by computing the gaugino loop in the presence of 
the TeV-brane gaugino mass M\, and then subtracting the value of the same diagram in the 
supersymmetric limit. The results are 

< = i6 j3+ i j2+ oH < i6 > 

< = i^+rH (17) 

< = M& + h Xi )> (18) 

2 W3_ 3 



™l = ^ 7^+™^ > (19) 



ml = -Uk), (20) 



2tt 2 V4 20" 

— (-' 

2tt 2 V5" 

and, if the Higgs doublets are effectively localized to the Planck brane, m\ u = m\ d = m~. Here 
X 3 , X 2 and X 1 are loop integrals defined as 



J a = J Q dq q 3 f Zt a[z = z = -; qj - f g>a \z = z' = -; q) 



M A =0 



(21) 



where a = 1, 2, 3 labels the gauge groups of the standard model. The functions f z , a (z, z'\ q) are 
defined in Appendix A (Eq. (57)). As shown there, the integrand in Eq. (21) has an exponential 
suppression ~ e~ 2q l T above T. This suppression at high momentum comes from the spatial 
separation between the matter fields on the Planck brane and the supersymmetry breaking on 
the TeV brane: to communicate the supersymmetry breaking to the MSSM scalars, the gaugino 
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in the loop must propagate from the Planck brane to the TeV brane, where its Majorana mass 
term is localized, and back. 

(iii) 321 KK states, XY gauge bosons, and XY gauginos 

We also calculate the tree-level masses of the 321 gauge-boson KK excitations, XY gauge bosons, 
and XY gauginos. For the 321 gauge bosons the equation determining the masses is 



J ° (t) + t m " Jl (a) _ Jo 



while for the XY gauge bosons we have 



(^) + (^) Y 



(22) 



J ( UhL\ J ( Hkt 



Yi (Y) Yo (^) ' 



(23) 



and for the XY gauginos, 



(24) 



In each case we will be interested only in the lightest massive mode. 

The results for the various spectra are shown in Figs. 2-4 for x ranging from .01 to 10. 
In these figures, we have chosen to normalize the particle masses, m, in units of 10 m g (i.e. 
the vertical axes represent 0.1m/mg), so that if the right-handed slepton masses are near their 
current experimental lower bound of around 100 GeV, the numbers labeling the vertical axes 
correspond roughly to masses in TeV units. In Fig. 2, the masses for the MSSM particles, XY 
gauge bosons and gauginos, and KK excitations of the 321 gauge bosons and gauginos are all 
shown. The value of T for a given x can be deduced from the fact that the XY and 321 KK 
gauge bosons have masses ~ 2.4 T; T ranges from ~ 100 TeV on the left hand side of the plot to 
~ 5 TeV on the right hand side. 3 From this plot one sees some overall features of the spectrum: 
the XY and 321 KK gauge bosons are quite heavy, with masses larger than 10 TeV even if one 
assumes that m~ e ~ 100 GeV. On the other hand the masses of the 321 gaugino KK excitations 
decrease much more rapidly as x is increased, and form pseudo-Dirac states with the lightest 
gauginos at large x, where the values of the masses plateau at m ~ T/4. Finally, the XY gaugino 
masses are lighter still, and continue to decrease with increasing x even after the 321 gaugino 
masses level off. In fact, we find that for large x the XY gaugino masses are given by [2] 



3 For these values of T, constraints from the precision electroweak measurements (T ^ 250 GeV [10, 14]) are 
completely negligible. 
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Figure 2: Masses of the MSSM scalars (dashed, with rriq, ma, and mg closely spaced and mj and 
mg below), MSSM gauginos (thick solid), XY gauginos (dot-dashed), 321 gaugino KK modes 
(dotted), and XY and 321 KK gauge bosons (thin solid, nearly degenerate and most massive). 
As explained in the text, we give the masses in units of 10 mg. 

where g 4 D = Qb/Vt^R is the "4D gauge coupling", which takes a value of 0(1). This fact 
significantly improves the discovery potential of the GUT particles at colliders. 

Figs. 3 and 4 focus on the low-lying masses. In Fig. 3 the masses of the MSSM scalars and 
gauginos are plotted. From the figure we immediately see that the ratios of gaugino masses to 
scalar masses become larger for larger values for x. This is because for large x (x ft 1) the 
scale of superparticle masses is close to the scale at which they are generated, so that the scalar 
masses are purely one-loop suppressed compared with the gaugino masses; on the other hand, 
for i C 1, the scalar masses are enhanced by a logarithm between the scale of superparticle 
mass generation and the gaugino masses, ln(l/x), and become close to the gaugino masses. 
Another interesting feature is that the ratios among the scalar masses are relatively insensitive 
to x for the range considered here, while the gaugino mass ratios change significantly as x does 
(their masses are less hierarchical for larger x, a feature that may actually be easier to see on 
the log-scale plot of Fig. 2). It turns out that the gaugino masses satisfy M a oc g\ for a; < 1 
and M a oc g a for x ^> 1 [11]. For scalar masses, we can see their rough scaling by studying 
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0.01 0.05 0.1 0.5 1 5 10 

x = M^/k 

Figure 3: Masses of the MSSM scalars (dashed) and gauginos (solid). 



the integrand of Eq. (21) evaluated at q ~ T, where the dominant contribution to the integral 
arises. This suggests that the scalar masses are approximately proportional to the square of the 
gauge couplings; in fact, we numerically find that the scalar masses scale roughly as the relevant 
Casimir times the relevant 4D gauge coupling squared {e.g. (4/3) (?| for the squark masses). 
Thus, the lightest among the superpartners of the standard model fields are the right-handed 
sleptons for the entire range of x considered. Including effects of the Yukawa couplings, the 
lightest one will be the right-handed stau. Since the lightest supersymmetric particle (LSP) of 
the model is the gravitino G 3 / 2 with mass of order T 2 /M P \ ~ .01 — .1 eV, it decays as f — > 
with a lifetime of order 8vrT 4 /m| ~ 1(T 18 - 1(T 14 sec. 

In Fig. 4 the masses of the XY gauginos and 321 gaugino KK modes are shown with the 
MSSM gaugino masses. We see that for x <: 0.4, the XY gauginos can be lighter than 2 TeV, 
well within the reach of the LHC. 4 For these larger values of x, the 321 gaugino KK modes also 
become light, and for much of the parameter space for which the XY gauginos can be discovered, 
these KK modes are also experimentally accessible. Therefore, the experimental signatures in 

4 Although the XY particles are strongly localized to the TeV brane, giving exponentially small couplings to 
matter (and ensuring sufficient proton stability), their production rates are still unsuppressed because they have 
order one couplings to the standard model gauge fields. 
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Figure 4: Masses of the MSSM gauginos (solid), XY gauginos (dot-dashed), and 321 gaugino 
KK modes (dotted). 

this parameter region are quite distinct: we would find two gauginos for each 321 group, whose 
masses are very close for large values for x, and one XY gaugino, which will be stable and 
seen as highly ionizing tracks at colliders [2]. 5 This raises the exciting possibility that both the 
underlying SU(5) gauge structure and N = 2 supersymmetric structure of the model will be 
revealed by discovering these gauginos at the LHC. 

It is interesting to point out that the three MSSM gaugino masses are determined by just two 
parameters, T and x, so that a single relation can be established among them. In particular, 
given the mass scale T, all the gaugino mass ratios are determined by a single mass ratio. 
This situation is depicted in Fig. 5, where we plot M 2 /Mi versus M 3 /Mi for the fixed value of 
rrie = 100 GeV. Moreover, the same two free parameters also determine the mass of the XY 
gaugino, M X y, so a measurement of, say, M 3 and M 1 would determine not only M 2 , but also 
M X y- The labeled points on the line in Fig. 5 are meant to give a sense of how M X y/M 1 varies 

5 In fact, if the theory preserves a certain parity of the bulk Lagrangian called GUT parity, the XY gaugino 
is absolutely stable. In the present case with matter strongly localized to the Planck brane, the XY gaugino is 
effectively stable for collider purposes even in the absence of GUT parity. 
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Figure 5: One gaugino mass ratio, M 2 /Mi , given as a function of another, M$/Mi_. Each point 
on the line specifies a different value of x — M x /k, and therefore a different value of M X y/M 1 . 
The labeled points correspond to Mxy/Mx = 0, 1, 3, 10, and 50. 

as the MSSM gaugino mass ratios change. These constraints among the gaugino masses make 
this part of the spectrum especially interesting from an experimental standpoint. In particular, 
one could predict the mass of the XY gaugino once the 321 gaugino masses were measured with 
sufficient accuracy. 6 

We finally discuss uncertainties. As was already mentioned in section 2.2, higher dimensional 
operators on the TeV brane introduce an uncertainty in the superparticle masses that we expect 
to be of order (/c/M*) 2 ^ 10%. We have used one-loop RGEs to evaluate the Planck-brane gauge 
kinetic terms at T, but errors coming from higher loops are small, especially because our results 
are not so sensitive to the relative size of the Planck-brane and bulk gauge couplings. Two-loop 
contributions to the sfermion masses could be sizable, especially for x <C 1: in this regime the 
logarithmically enhanced higher-loop effects become larger. For x ^ .01, however, we expect 
that those are still O(10%) effects. Based on these arguments we expect that our results are 
robust, at worst, at the 20 — 30% level. 

6 It is worth mentioning that the gaugino mass relations shown in Fig. 5 (and in Fig. 4 up to an overall mass 
scale) are insensitive to certain model details. In particular, these results apply even if matter propagates in the 
bulk, in which case squark and slepton masses can arise at tree level, allowing much lower values for T. 
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4 Electroweak Symmetry Breaking 



The naturalness of the electroweak symmetry breaking sector is an important issue for any ex- 
tension of the standard model, and the original motivation for low-energy supersymmetry. The 
fact that we have not seen physics beyond the standard model typically requires some degree 
of fine-tuning among parameters in this sector. In supersymmetric theories, the experimental 
results that are most relevant for naturalness considerations are (i) the non-discovery of super- 
particles, and in particular gauginos, which can easily be lighter than the other superparticles 
because their masses are protected by an R symmetry, and (ii) the non-discovery of the lightest 
Higgs boson, which is predicted to be light in supersymmetric extensions of the standard model. 
In the MSSM, the physical mass of the lightest Higgs boson can be as large as ~ 130 GeV, 
but for the Higgs boson to have evaded detection at LEP II requires somewhat large top squark 
masses. These in turn generate a large negative contribution to the soft mass-squared parameter 
for the Higgs field, introducing some degree of fine-tuning in the Higgs potential. 

In our theory all superparticles are heavy enough to evade experimental bounds provided 
that the right-handed slepton mass is larger than about 100 GeV. This is clear from Figs. 2 
and 3: the gaugino masses are significantly heavier than the scalars, especially for larger x, 
because they acquire tree-level masses. The right-handed stau will be somewhat lighter than 
the other right-handed sleptons due to Yukawa-induced radiative corrections (not included in 
our analysis), but the experimental bound on the stau mass is also somewhat less stringent, so 
we keep 100 GeV as a representative value for the lower bound. 

Now we can ask the following question: for a given value of x, if we take T large enough 
to evade the bound on the right-handed slepton mass, how heavy is the lightest Higgs boson? 
Keeping all other parameters fixed, the Higgs mass decreases as x increases because the top 
squark masses decrease (see Fig. 3). The left-right mixing in the top squark mass matrices is 
dominated for moderately large tan (3 by a loop-generated A t term, 



Here we have included only the SU(3)c loop, yt is the top Yukawa coupling, and h z $ is the 
function defined in Appendix A (for SU(3)c)- We find that the mixing induced is quite small. 
A no-mixing scenario requires a somewhat larger overall mass scale for superparticles, especially 
for x <; 1, to obtain a sufficiently large Higgs boson mass: we estimate mg <; 200 GeV is required 
for x 1. Even with these large values of superparticle masses, however, the fine-tuning in our 
theory is not as severe as one might naively imagine based on the squark masses alone. This is 
because the scale of the superpartner masses is close to the the scale where they are generated, 
especially for large values of x, and thus even if the top squark is somewhat heavy, it will not give 




(26) 
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too large a negative contribution to the Higgs mass-squared parameter. This is an interesting 
point: the correction to the physical Higgs-boson mass arises below the stop mass, and goes 
like \n(mi/m t ), while the correction to the Higgs mass-squared parameter arises above the stop 
mass, and goes like ln(M mess /m^), where M mess is the scale where soft masses are generated. 

An alternative way of obtaining a large Higgs boson mass is to introduce a singlet superfield 
S on the Planck brane with the superpotential interactions 25(i/)(\SHdHd + kS 3 ). This setup 
is also motivated as a way to naturally induce the supersymmetric mass term (// term) for 
the Higgs doublets through the vacuum expectation value for the S field. 7 In this case the 
Higgs boson mass receives additional contribution at tree level thorough the coupling A. The 
size of this contribution depends on the value of A, whose upper bound is set by Landau pole 
considerations. An interesting point is that in our theory the 321 gauge couplings become strong 
at high energies so that the bound on A is significantly weaker. (This fact can be understood 
more easily in terms of the 4D dual picture of the theory [2].) This allows us to have the weak 
scale value of A as large as A ~ 0.8, and we obtain a large enough Higgs boson mass for tan (3 ^ 5 
even with mg ~ 100 GeV. 

How finely tuned is electroweak symmetry breaking in our theory? A precise discussion 
of naturalness requires the values of the Higgs mass-squared parameters, and these depend on 
the details of the Higgs sector. Here we simply estimate their sizes for Planck-brane localized 
Higgs doublets, in which case the soft masses vanish at tree level. The dominant radiative 
corrections come from the one-loop SU (2) l gauge contribution and the two-loop top- Yukawa 
S'L / "(3)(7-gauge contribution. The former is calculated in the previous section (given by Eq. (19) 
with m\ u \ a2 = m 2 hd \ a2 — m~), which gives m\ u \ a2 ~ 4m? for large x. The latter we expect is 
similar to the flat space case for x 1, with the gluino masses in the two theories identified. 
Referring to [15], we obtain mfc \ a3at ~ — 5m|, although the precise value is quite sensitive to the 
choice of the renormalization scale etc., and the exact number we quote is not really trustworthy. 
Nevertheless, it is reasonable to expect m\ = m\ \ a2 + Tn\ \ asat — ~ cm % where c is a factor of 
a few. The amount of fine-tuning is given roughly by |m^/m/ lu | 2 /2 for moderately large tan/3, 
and this is of order 10% if we take m g at its present experimental bound. This is relatively 
mild. The bound on the physical Higgs boson mass may push up the superparticle mass scale 
higher, and for smaller x, the fine-tuning might change due to ln(l/x) enhancements in the loop- 
induced scalar masses. However, we still expect that the fine-tuning is not very severe. Although 
this estimate is based on top-derived radiative electroweak symmetry breaking, it illustrates the 
required fine-tuning in our theory even for more general Higgs sectors. 

7 With these superpotential interactions, the U(1)r symmetry discussed in section 2.1 is explicitly broken to 
the Z4.R subgroup, under which S carries a charge of +2. This Z^,r symmetry, however, is still sufficient to 
forbid unwanted dimension four and five proton decay operators and a large supersymmetric mass term for the 
Higgs doublets. 
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We finally comment on the possibility of a very simple and constrained scenario. Suppose that 
both matter and Higgs fields are localized to the Planck brane (either exactly or approximately, 
by bulk hypermultiplet masses), and a ji term of order the weak scale is generated on the 
Planck brane. Then, assuming top- Yukawa driven radiative electroweak symmetry breaking, the 
entire superpartner spectrum and the parameters of the Higgs potential can all be calculated in 
terms of only three parameters, x, T and /x, one of which is fixed by the observed value of the 
Higgs expectation value v. In this setup, not only supersymmetric flavor problem, but also the 
super symmetric CP problem is solved, because the A terms and the B term are all generated 
radiatively through gaugino loops, and are thus all real in the basis where the gaugino masses 
are real. The sign of \x is determined to be negative in the standard phase convention. A detailed 
study of this scenario will be interesting. 

5 Conclusions 

Warped supersymmetric grand unification has some remarkable features. Most notably, it pre- 
serves MSSM-like gauge coupling unification even though it predicts that the exotic particles 
associated with grand unification - XY gauge bosons, for example - appear near the TeV scale. 
These light exotics would normally induce baryon-number violating processes at disastrous rates, 
but here proton decay is naturally suppressed simply by localizing matter to the Planck brane. 

This framework also accommodates simple, constrained possibilities for supersymmetry break- 
ing. In this paper we considered a setup with unsuppressed, or only mildly suppressed, super- 
symmetry breaking localized to the S , t/(5)-preserving TeV brane. In such a scenario, the su- 
perpartners of the standard model fields naturally acquire TeV-scale masses. Taking the quark 
and lepton superfields to be localized to the Planck brane, advantageous for suppressing proton 
decay and enforcing flavor universality in the squark and slepton masses, one is led to a highly 
predictive model. In this case, the masses of the MSSM squarks, sleptons, and gauginos, along 
with those of the XY gauge bosons and gauginos and KK excitations of MSSM gauge particles, 
are all calculable in terms of two parameters: T, the scale of the infrared brane, and x = M\/k, 
the ratio of the supersymmetry-breaking gaugino mass on the TeV brane to the curvature scale. 
Calculating these masses was our main purpose. 

In this setup the gauginos acquire masses at tree level and are heavier than the squarks and 
sleptons, which only acquire masses at loop level, and the lightest MSSM superparticles are the 
right-handed sleptons (including Yukawa effects the lightest one is the right-handed stau, which 
decays into the LSP gravitino). Requiring that the right-handed sleptons are heavy enough to 
evade detection at colliders sets a lower bound on T for a given value of x. From this lower 
bound we infer that the XY gauge bosons and KK excitations of the standard model gauge 
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bosons are too heavy to be detected at the LHC, with masses larger than 10 TeV regardless of 
how large x becomes. On the other hand, supersymmetry breaking has a dramatic impact on 
the masses of the supersymmetric partners of these gauge bosons, pushing some of the gaugino 
states to be considerably lighter. In fact, for x <^ 0.4, we find that the XY gauginos can have 
masses below 2 TeV, and for these larger values of x, the KK excitations of the MSSM gauginos 
are also relatively light, and approach degeneracy with the MSSM gauginos themselves as x is 
increased. The fact that the masses of these particles can be light enough to be discovered at 
the LHC - which would reveal both the underlying unified gauge symmetry and the enhanced 
N = 2 supersymmetry of the theory - is the most important result of this paper. 

Our results for the spectra of GUT particles and superparticles are independent of the Higgs 
sector, but we also briefly considered electroweak symmetry breaking in the particularly con- 
strained setup where the Higgs doublets are localized to the Planck brane. In the minimal 
scenario, the experimental bound on the Higgs boson mass may require the right-handed slep- 
ton mass to be somewhat larger than its experimental lower bound of 100 GeV. On the other 
hand, a superpotential coupling between the Higgs doublets and a Planck-brane localized singlet 
S, motivated independently as a means for generating a weak-scale \i term, allows the Higgs mass 
to be raised above its lower bound quite easily. The effect on the Higgs mass can be stronger 
than in the conventional NMSSM because our theory has relatively strong gauge interactions at 
high energies, which drive the SHe>He> coupling away from its Landau pole. It will be interesting 
to explore electroweak symmetry breaking in this model in greater detail. 
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Appendix A 



A.l Gaugino propagators 



In this appendix we derive the propagators for the 321 gauginos in the presence of general brane 
kinetic terms and a TeV-brane localized Majorana mass. The free action for the gauginos is 
given by 



S = 
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where A and A' are the two-component gaugino fields that are contained in 4D superfields V and 
S, respectively. Here, we have suppressed the index a running over SU(3)c, SU(2) L and U(l)y 
for the simplicity of the notation, but it should be understood that the Planck-brane localized 
kinetic term, g , takes different values for SU(3) C , SU(2) L and U(1) Y - go 9o,a- The other 
parameters, gs, 9-k and M x , are SU(5) symmetric and have universal values for SU(3)c, SU(2) L 
and U{1) Y . 

It is useful to define the rescaled gaugino fields A = e~ 2ky X and A' = e~ 2ky X' . Then, in terms 
of the Fourier transformed fields X(p,y) = j d A xX(x 1 y)e ipx and X'(p,y) = J d A xX'(x,y)e ipx , the 
above action is written as 
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The content of this action can be divided into two parts. First, by examining the region near 
y = and nR in the equations of motion derived from Eq. (28), we find the following conditions 
imposed on the fields: 
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together with the equations obtained by interchanging A(p) <-> A^(— p), A'(p) <-> A't(— p) and 
j n the above equations; here e — > 0. These conditions provide boundary conditions for 
the corresponding propagators because the propagators are given by the time-ordered products 
of the two fields evaluated on the vacuum: G w > = (0\T{(f(p'}\0). Second, the bulk piece (the 
terms that do not involve delta functions) of Eq. (28) is written as 



Sbui k = lj4^C d y(^(P) A'(-p)|A(-p) A't(p)) 
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(33) 



This piece dictates the form of the propagators in the bulk. Defining the 2x2 matrix appearing 
in Eq. (33) as M bulk , the propagators 
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are given as a solution of 
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in the bulk, where 1 is the 4x4 unit matrix. Note that G represents propagators for the 
rescaled fields, A and A', which are related to propagators G for the unrescaled fields, A and A', 

as G = e 2k( -y + y'^G. 

Now, let us derive the bulk propagator by solving Eq. (35). Parameterizing the matrix G of 
Eq. (34) as 



G = 



i^Pv f 


ie~ ky (d y - §; 
ia^Pf, f 


r 


i h 

^e^( dy 


\)h 


^e~ k y(dy 

ih' 


\)h' 


i h 

^e- k y(-d y - 1) h 


ih' 


\)h> 


i^Pn f 
ie- ky (-d y - \ 


)/ 


ie- k y{d y - f 

ia^Pf, f 


)/' 



(36) 

where /, /', h and h' are the functions of y, y' and p, we find that the functions /, /', h and h' 
obey the following equations: 



1 

7b 



pV» + e -ky {d 2 _ kdy _ k 2 } f ^ y '. p) = S (y - y'), 



(37) 
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9l 



p ^y + e -kv {d 2 _ kdy + , 2) ^^jp) = % _ ^ 



p2e*» + e~ k y(d 2 y - fcfy - -A; 2 ) h(y,y';p) = 0, 



4 

Changing the variable from y to z = e ky /k, the equations become 

f z (z,z';p) = 5(z-z'), 

h z (z,z';p) = 0, 
K(z,z';p) = 0, 



1 


~P 2 + dl- 


3 1 " 

4^2 


i 


p 2 + d 2 z + 


1 1 ■ 

4^2 




p 2 + ^ 2 - 


3 1 ■ 

li 2 ". 




p 2 + d 2 z + 


1 1 ■ 

li 2 ". 



(38) 
(39) 
(40) 

(41) 

(42) 

(43) 
(44) 



where the functions with the subscript z represent the functions obtained by changing variables 
from y and y' to z and z'\ f z (z,z';p) = f(\n( y kz)/k,\n(kz ; )/k;p) and so on. Considering the 
region z ^ z', we find that the solutions to Eqs. (41 - 44) are given as 



fz(z,z';p) 
f z (z,z';p) 
h z (z,z';p) 
h' z (z,z';p) 



= sfr Wiz') h(\p\z) + b > (z') KMz)), 



z\^Wl (\p\z) + K (\p\z)), 



^(a^z') h(\p\z) + P>(z') Mlplz)) , 
sTzU'^z') h{\p\z) + (3'^z') K (\p\z)) , 



(45) 
(46) 
(47) 
(48) 



for z > z', where I n (x) and K n {x) are the modified Bessel functions of order n and \p\ = \J — p' 2 . 
For z < z' the functions a > , b > , a>, &>, a>, /3>, a> and /3> take different forms, which we denote 
as a<, 6<, a<, &<, a<, /3<, a< and respectively. 

The propagators must obey the conditions of Eqs. (29 - 32). This allows us to solve the 
functions 6>, &>, /3> and /3> in terms of a>, a>, a> and a>; and similarly, 6<, &<, /3< and /?< 
in terms of a<, a' a< and a<. Further constraints on the remaining functions come from the 
continuity of the propagators at z — z'\ 



fz(z,z';p)\ z=z , +e 

f'z( Z , Z ';P)\z=z'+e 

h z (z,z';p)\ z=z > +e 
h' z {z,z';p)\ z=zl+e 



fz(z,z';p)\ z=z/ _ e , 

f'z( Z , Z ';P)\z=z>-e, 

h z (z,z';p)\ z=z >- e , 

K( Z i Z ''iP)\z=z'-e, 



(49) 
(50) 
(51) 
(52) 
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and the junction conditions following from Eqs. (41 - 44) at z — z'\ 

d z f z (z,z';p)\ z=z , +£ - d z f z (z,z';p)\ z=z ,- e = g 2 B , (53) 

dzf' z (z,z';p)\ z=z , +e - d z f' z (z,z;p)\ z=z/ - e = g 2 B , (54) 

d z h z (z,z';p)\ z=z , +e - d z h z (z,z';p)\ z=z ,_ e = 0, (55) 

d z ti z (z,z';p)\ z=z , +t - d z h' z (z,z';p)\ z=z/ - e = 0. (56) 

These equations completely determine the functions a>, a^, a>, «>, a<, a^, a< and a' K . 

After some algebra, we finally find that the gaugino propagators defined in Eq. (34) are given 
by Eq. (36) with 

U^z'-p) = { f^^ B2 (*i(bl*<) + CK^pM) 

x ({C - A){h(\p\ z> ) + AKMz*)} - B'KMz^ , (57) 

f,Mp) = (*>(bl*<) - ck {\p\zS) 

x ((C - A){/ (b|2>) - Atf (bl*>)} + B 2 K {\p\ Z> )^ , (58) 

= - { c B -^AfT& { IMz<) + ck 1 {\p\z<))b[i 1 {\p\z > ) + ck 1 (\ p \z > ) ) ), (59) 
= { c-^+b^ { IMz<) ~ c^o(bk<))g(/o(bl^) - c/r (bl*>)), (60) 

where |p| = ^f—p 1 and z< (z>) is the lesser (greater) of 2 and z'; the functions f z , f' z , h z and h' z 
are related to /, /', h and ft/ by f z (z,z';p) = f(y,y';p) etc. with z = e ky /k. The coefficients A, 
B and C are given by 

, XjX K — YjY K XiYk + XkYi „ Zj 

A = K n , -D = 7i n , = . 61 



Here, 



^ = ^o(¥) - l ^ = M A ^(M), | = ^o(M) + |^(¥), 

(62) 

where T = fce _7rfcil . 



A. 2 Limiting behaviors 

We here consider various limits of the gaugino propagators. We first consider the "4D limit" , 
in which the compactification scale is sent to infinity and the theory reduces to the 4D MSSM. 
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This limit is obtained by taking \p\/T — > keeping T//c and |p|/M A fixed. The function f z then 
becomes 

Uz ' z ' ;p) ~-^w + (iM^r (63) 

where g 4 D is the 4D gauge coupling given by l/<?f D = TtR/g 2 B + 1/(?q + l/<7^. The MSSM 
gaugino propagator, G^(p), is obtained by multiplying the propagator in this "4D limit", 
G AA t (y, y'; p) = icr tJ, p lx f(y, y'; p), by the MSSM gaugino wavefunction in X(x, y), e 3ky / 2 , and setting 
y — y' — 0. Considering the 4D gaugino mass is given by Mx,4D = 9\t>M\T /k (see Appendix B), 
this is given by 

G%(V)=9^-^^, (64) 

P 1V1 X,4D 

reproducing the 4D gaugino propagator. 

We next consider the propagator with the external points both on the Planck brane in 
the limit where the momentum scale \p\ is much larger than the IR scale T: \p\ 3> T. We 
are interested in the difference between the gaugino propagator in the presence and absence of 
supersymmetry breaking. For \p\ ^> T, the relevant quantity f(p) = f z (z — z' — l/k;p) — f z (z = 
z' = l/k;p)\ Mx =o is given by 

/ f 1 (l£l)JC (lEi)+J (lEl)jf 1 (l£i) 

which shows that f(p) is exponentially suppressed for p> T. This ensures the UV insensitivity 
for the scalar masses computed in section 3 (see Eqs. (16-21)). 

Appendix B 

In this appendix we present a formula for the 321 gaugino masses, applicable for any value of 
x = M\/k, in the presence of general brane-localized kinetic terms. We start with the action 
given in Eq. (27). This action gives the following equations of motion in the bulk: 

p ky ~ 1 k * 

—ia^X + —(-d y + -)A't = 0, (66) 
9b 9b * 

p ky - 1 k - 

— «7*0„A* + -(d y + -)A = 0, (67) 
9b 9b z 
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where we have presented the equations in terms of the rescaled gaugino fields, A = e 2ky \ and 
A = e~ 2ky \'. Looking for solutions of the form 



X(x,y) = J2\ n (x)ti(y), X'(x,y) = £ A n (x) /*'(</), 



(68) 



the bulk equations of motion, Eqs. (66, 67), lead to the following differential equation for /, 

3A; 



A. 



d 2 X ~ kd vfn ~ -rti + rny k yf* = 0. 



(69) 



Here, m n is the 4D masses and we have used the 4D relation ia^d^\ n = m„Aj,. The solution of 
this equation is given by 



f x n {y) 



Ji 



k 



-Mi 



+ b x Y 1 



~k 



(70) 



where N n and 6a are coefficients that do not depend on y. 

The boundary conditions for f£ and aty = and tcR are given by examining the equations 
of motion, Eqs. (66, 67): 



eX' 



, 9b r\ 
y=t 9o 

+ (dy + \)f. 



A 



0, 



fA' 



y=e 

y=7rR—e 
nkR f X' 



o, 



n 2 „nkR 

+ 9 -^m n r n 

9l 



y=e 
A 



y=nR 

+ (dy + \)fn 
y=nR—e £ 



- glMxfi 



y=irR 



o, 



y=irR—e 



= o, 



(71) 
(72) 
(73) 
(74) 



where e — > 0. Eliminating gives boundary conditions for at y = and 7ri2, each of which 
determines the coefficient b\. Equating b\ obtained from the boundary conditions at y = with 
that from y = irR, we obtain the equation that determines the 321 gaugino masses: 



JI 

lo (^) + fm^ (^) F - f |m n F! + g 2 B M x Y 1 

where T = fce - ^. 



m n Ji 



T 



+ ^|M A J 1 (^) 



(75) 
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